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QCD dynamics at small quark and gluon momentum frac- 
tions or large total energy, which plays a major role for HERA, 
the Tevatron, RHIC and LHC physics, is still poorly under- 
stood. For one of the simplest processes, namely bb pro- 
duction, next-to-leading-order perturbation theory fails. We 
show that the combination of two recently developed theoret- 
ical concepts, the fcx-factorization and the next-to-leading- 
logarithmic-approximation BFKL vertex, gives perfect agree- 
ment with data. One can therefore hope that these concepts 
provide a valuable foundation for the description of other high 
energy processes. 

Existing QCD calculations describe many high energy 
observables which involve partonic transverse momen- 
tum rather poorly. This is also true for the theoreti- 
cally especially clean case of bb production, which was 
investigated experimentally at Fermilab |Q . Since central 
quark- antiquark production at y/s = 1.8 TeV is sensitive 
to very small gluon momentum fraction x « 10~ 2 — 10~ 4 , 
one probes the gluon content of the nucleon at small x, 
which is a central issue of current research. We recon- 
sider this process and combine as essential new ingredi- 
ents the fc t -factorization scheme with the next-to-leading- 
logarithmic-approximation (NLLA) BFKL production 
vertex derived in H. The fct-factorization approach 
for the description of high energy processes |^ |(| differs 
strongly from the conventional NLO collinear approxima- 
tion (e.g. 0) because it takes the non- vanishing trans- 
verse momenta of the scattering partons into account. 
The usual gluon densities are replaced by unintegrated 
gluon distributions which depend on the transverse mo- 
mentum fc t . These together with the fc t -factorization 
form a basis for a general calculation scheme for high 
energy (i.e. small x). The standard collinear approxi- 
mation has the advantage of being closely related to the 
operator product expansion. It is, however, only justified 
for the processes dominated by x = 0(1). In application 
to processes governed by small x the k t -factorization ap- 
proach has the advantage that its approximations corre- 
spond to the dominant kinematics. Essential small x con- 
tributions are included in the Born approximation which 
in the collinear approach are accounted for in higher or- 
ders only. This is well known from the case of structure 
functions where the DGLAP evolution is appropriate for 
x = 0(1) and the BFKL evolution for small x. 

While the fc f -factorization formalism is very attractive 
theoretically, its phenomenological usefullness has been 



mostly tested in the case of the structure function F% 
§,|o|. The NLLA BFKL vertices are just the ones needed 
to treat semi-hard central production at collider energies 
in this approach. 

In our calculation we use one particular element of the 
NLLA BFKL formalism |QJ|, namely the effective ver- 
tex for quark-antiquark production. Thus our calcula- 
tion can be seen as a first phenomenological application 
of this vertex which decides whether the NLLA BFKL 
formalism can be hoped to converge. 

One special aspect of the reaction we investigate is the 
possible loss of gauge invariance when a qq production 
vertex is incorporated into an amplitude with off-shell 
gluons. In the BFKL approach, however, gauge invari- 
ance is ensured automatically by the use of the just men- 
tioned NLL effective vertex which is valid in quasi multi 
Regge kinematics (QMRK), i.e. when the q and q have 
similar rapidities and form a cluster (in contrast to LLA, 
where the particles are produced with a large rapidity 
gap). 




FIG. 1. The basic diagram 

We begin with the following definition for the light cone 
coordinates and the momenta of the scattering hadrons 
in the cm. frame 

fc+ = fc° + fc 3 , k~ = k° - fc 3 , k x = (0, fc 1 , fc 2 , 0) = (0, k, 0). 
P? = P 2 ~ = y/i, Pi = Pt = 0, Pi± = P2± = 0. 

The Mandelstam variable s is as usual the cm. energy 
squared. As defined in fig. |l|, qi and q 2 are the momenta 
of the gluons and the on-shell quark and antiquark have 
momentum k\ respectively ki- In the high energy (large 
s) regime we have 
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The longitudinal momentum fractions of the gluons are 

x i = q.i l p i, x 2 = -<h l p 2 ■ 

The cross section for heavy quark pair production in 
the fct-factorization approach is then given by |]H 
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The factor (N 2 — l) 2 reflects the projection on color sin- 
glet, where N is the number of colors. The hard ampli- 
tude ip C2Cl (x%,X2, qii_, q 2 ±, kx,k2) is calculable in pertur- 
bation theory, whereas the unintegrated gluon distribu- 
tion T(x, q±) has to be measured or modelled. We choose 
the argument [x 2 of the strong coupling constant as{^ 2 ) 
in the hard amplitude ip C2Cl to be equal to q 2 = — q 2 ± 
respectively q| = 0. 
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FIG. 2. The effective vertex 

We generalize the results on the qq production vertex 
presented in Q for massless QCD in an obvious way in 
order to take the masses m of the produced quarks into 
account. The resulting vertex ^ C2C1 is given by a sum of 
two terms 

* C2C1 = -g 2 (t c H C2 b(k 1; k 2 ) - t C2 t c > b T (k 2 , kx)) , (2) 

where t c are the colour group generators in the funda- 
mental representation. The connection between ip C2Cl in 
eq. ([!]) and 'F C2Cl in eq. (^) is given by 

ip C2C i = u(fci)^ C2Cl v(fc 2 ), 

with the on-shell quark and antiquark spinors u(fc) and 
v(fc). The expression for b(ki, k 2 ) is a sum of two terms 

w, , x -Aii--jki±-m + lpT+~ f) (q 2 ,q l ) 

b(k 1 ,k 2 )=j - — - 7 + — , 

{qx — kiY — m^ {k\ + k 2 ) A 

(3) 

The first term on the r.h.s. of eq. (^) describes the pro- 
duction of a qq pair by means of usual vertices (see fig. |^), 



the second term involves the light-cone projection of the 
effective vertex T + ~ /3 (q 2 ,qi), which describes the tran- 
sition of two t-channel gluons (reggeons) with momenta 
qi and q 2 to a gluon with momentum k\ + k 2 
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with ti/ 2 = q\j 2 - This effective vertex differs from the 
usual triple-gluon vertex by the appearence of the last 
two terms. They are related to Feynman diagrams in 
which the qq pair is not produced by the ^-channel gluons 
but in other ways. These two last terms in eq. (|J) are also 
required by gauge invariance, T + ~ ^(q 2 , <7i)(qi —<l2)f3 = 0. 
Another consequence of gauge invariance is the vanish- 
ing of the matrix element of the effective vertex ^ C2 Cl 
between on-mass-shell quark and antiquark states in the 
limit of small q±± or q 2 ± 

u(fci)^ C2 Cl v(fc 2 ) -> for q 1± or q 2A _ -> 0. 

The function b T (k 2 ,ki) is very similar to (||) 



b T (k 2 ,k 1 ) =Y 



Ai-L- M± + m _ ypT + /3 (g 2 ,9i) 
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The unintegrated gluon distribution is related to the 
standard gluon distribution by 



,q 2 ) = j (°°^e( q 2 -k 2 )^(x,k). 



xg(x 



Taking the derivative of this expression makes it obvious 
that T{x, k) includes the evolution of xg(x, q 2 ), which is 
given by the BFKL and/or DGLAP equation. Since the 
unintegrated gluon distribution is not known at small k, 
we write this equation as 



xg(x,q 2 ) = xg{x,ql) + 



dkj 
k 2 



9(q 2 -k 2 )^(x,k). (5) 



This formula has been repeatedly used |IJ,|l3|,f|,[l4| and 
introduces the a priori unknown initial scale go and the 
initial gluon distribution xg(x,q 2 ). Following one 
may neglect the hard cross section dependence on q in 
the soft region |q| < go, so that 

~ 1 ~TTj2 — Tvf f ~2~ = z{qix,q2X) ->■ 

5(g 1± ,g 2± )e( q 2 -g o 2 )0(q2-9o) 
+ 5(g 1± ,0)e(q 2 -g 2 )6(g 2 ~q 2 ) 
+ 5(0,g 2 i)e(q 2 -g 2 )e(g 2 -q 2 ) 
+ 5(0,0)e(g 2 -q 2 )e(g 2 -q 2 ). (6) 

Note that the very existence of the finite limit q± — > 
follows from the decrease of the production amplitude 
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due to gauge invariance. Substituting this formula in ([!]) 
using eq. (^) one may easily perform the integration over 
q±. As a result, 5(0, 0) produces the standard expression 
of collinear factorization (ref. [jl0],[l4]] ) , while S(qi±,0), 
5(0, <72_l) correspond to the asymmetric configurations, 
where one of the gluons is described by the unintegrated 
distribution and the other by the integrated one. Here 
it is important to notice that when we insert (|5|), (|^) in 
(|jj) the coupling constant a s in the term proportional to 
xg(x, q 2 ) is taken to be a s {qo). 
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FIG. 3. The calculated b cross section in comparison to 
experimental data from the DO Collaboration and the NLO 
QCD result with MRSR2 structure functions and nib = 4.75 
GeV from [15] 

In all our numerical calculations we used for the uninte- 
grated gluon distribution T{x, k) the code by Kwiecinski, 
Martin and Stasto |l(|, because they use a combination 
of DGLAP and BFKL equations which governs simul- 
taneously the evolution in Q 2 and x. They obtain an 
excellent description of ^(x, Q 2 ) in a very large x-Q 2 - 
window. According to our knowledge this is the only 
unintegrated gluon distribution which has given such a 
satisfactory result, which justifies our choice. As in the 
case of the usual gluon distribution function one has to 
choose an initial scale and an initial distribution function 
which in the case of [nO| are given by 



ql = 1 GeV 2 , xg(x, q 2 ) = 1.57(1 - x) 



2.5 



(7) 



We use these values, which are fixed by the fit to 
F 2 (a;,Q 2 ), in our calculation. 

We consider the production of 6fe-pairs. For the com- 
putation we use eqs. (|l|) and (||) with the unintegrated 
gluon distribution from JlCfl and the corresponding val- 
ues The rapidities and the transverse masses of the 
produced quark and antiquark are defined by 



2/1/2 = \Mt=^)> m l/2± 
Z /C l/2 



k 2 

ft l/2_l_- 



The Bj or ken- variables of the gluons can then be written 
as 

Xi = -^{m 1± e Vl +m 2 ±e y2 ), 

x 2 = —(m x ±e~ Vl +m 2 ±e~ y2 ). 
\ s 



In fig. H we show our results for inclusive b production, 
together with experimental results measured by the DO 
Collaboration @ (see Table II) in^i =1.8 TeV pp col- 
lisions. We obtain this cross section by integrating out 
all antibottom variables in cq. (Til). 
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FIG. 4. The result for the semi differential 66 cross section 
at fcix m i n = 6.5 GeV, compared to CDF data and the NLO 
QCD result with MRSD0 structure functions and mt = 4.75 
GeV from [1] 

The variable fci±min is the lower integration cut on the 
transverse momentum of the produced b quark. To get 
an indication of the theoretical uncertainties apart from 
higher order contributions which are not available at the 
moment we proceed in a similar way as the authors of ref. 
jl| and present our calculations for three different choices 
of Aqqjj) and the bottom quark mass 
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A (5) = 180 MeV, m b = 4.5 GeV, 
A (5) = 150 MeV, m b = 4.7 GeV, 
A (5) = 100 MeV, m b = 4.9 GeV, 



Our result is in very good quantitative agreement with 
data over the whole range of fcij_ m i n . The corresponding- 
central QCD NLO calculation has a similar shape, but 
is about a factor of 2 — 3 smaller than our central result 
(see for example fig. 11 in [p§). 
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5. The result for the semi differential bb cross section 
= 8.75 GeV, compared to CDF data 



We now turn to bb correlations in y/s = 1.8 TeV pp 
collisions, which have been measured by the CDF collab- 
oration at Fermilab . The correlations of the quark and 
antiquark give an insight into the dynamics of the pro- 
duction mechanism and are important in order to study 
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the limits of the collinear (k±± = — &2_l) LO QCD approx- 
imation. We present a comparison between our results 
and the experimental data in fig. [| and fig. 0. The data 
points and uncertainties were taken from jj],[17| . We find 
good agreement with experiment for both fcixmin = 6.5 
GeV (fig. |) and fci_L min = 8.75 GeV (fig. |). In this 
case QCD NLO calculations underestimate the measured 
cross section roughly by factor of 3 (compare with fig. 6 

0,1 



-e- 

"C 1E-4 



Sqrt(s) = 16TeV 




|y,J<i ^ 


— central 
high 
low 




Sqrt(s) = 1.8TeV 



20 40 



80 100 120 140 160 180 

(|) [Degrees] 

FIG. 6. 4> distribution of bb hadroproduction 



ization gives them only at NLO and beyond. In contrast, 
the available NLO caculations jl8| are not in agreement 
with the Tevatron data we compare with [p] JT^Jl^] . 

Our results show that at least those features of the 
effective qq vertex which we tested provide a substan- 
tial improvement with respect to the standard collinear 
treatment. 

If further tests of other observables should be equally 
successfull, the NLL BFKL vertices will also allow for 
a much improved description of many processes to be 
studied at RHIC and LHC. 

We thank J. Kwiecinski, A. Martin and A. Stasto for 
supplying us with their code for the calculation of the 
unintegrated gluon distribution. 

L. S. would like to acknowledge the support by the 
DFG. 



An interesting parameter concerning the correlation is 
the opening angle <j) between the momentum vectors of 
the produced quarks in the plane transverse to the beam 
axis. Our predictions for the corresponding differential 
cross sections at Fermilab and LHC energies are shown 
in fig. H As expected we find a peak at <f> = 180° which 
shows the dominance of the collinear part. 
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FIG. 7. Rapidity distributions of bb hadroproduction 

Additionally we present our predictions for rapidity 
distributions of the b for the rapidity of the b being and 
y/~s = 1.8 TeV respectively ^/s — 16 TeV in fig. ^ Our 
cross section for = 1.8 TeV at j/2 = is about a factor 
of 3 larger than the corresponding QCD NLO result from 

Let us conclude. We have studied quark-antiquark 
hadroproduction within the /c t -factorization approach us- 
ing an unintegrated gluon distribution and a specific ef- 
fective BFKL vertex for qq production. We found very 
good agreement with experiment for both single b pro- 
duction and bb correlations at yfs = 1.8 TeV. Our ap- 
proach leads to nontrivial bb correlations already at LO 
perturbation theory, whereas traditional collinear factor- 
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